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X-ray spectral variability in the ultraluminous X-ray 
source Holmberg IX X-1 



0^ 
O 
O 

o 

Q 
in 



Oh; 

o 



> 

o 

a^ 

(N 

On 
O 



K. Vierdayanti^*, C. Done^, T. P. Roberts^, and S. Mineshige^ 

^Department of Astronomy, Kyoto University, Kitashirakawa Oiwake-cho, Kyoto 6068502, Japan 
^Department of Physics, University of Durham, South Road, Durham DHl 3LE, United Kingdom 



Accepted 2009 December 00. Received 2009 December 00; in original form 2009 October 00 



ABSTRACT 

We use XMM-Newton and Swift data to study spectral variability in the ultraluminous 
X-ray source (ULX), Holmberg IX X-1. The source luminosity varies by a factor 3-4, 
giving rise to corresponding spectral changes which are significant, but subtle, and not 
well tracked by a simple hardness ratio. Instead, we co-add the Swift data in intensity 
bins and do full spectral fitting with disc plus thermal Comptonisation models. All the 
data are well-fitted by a low temperature, optically thick Comptonising corona, and 
the variability can be roughly characterised by decreasing temperature and increasing 
optical depth as the source becomes brighter, as expected if the corona is becoming 
progressively mass loaded by material blown off the super-Eddington inner disc. This 
variability behaviour is seen in other ULX which have similar spectra, but is opposite 
to the trend seen in ULX with much softer spectra. This supports the idea that there 
are two distinct physical regimes in ULXs, where the spectra go from being dominated 
by a disc-corona to being dominated by a wind. 

Key words: accretion, accretion disc - black hole physics: stars: individual: Holmberg 
IX X-1 - X-rays: binaries. 



1 INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are X-ray sources that 
are too bright (Lx > 10''® ergs s~^) to be powered by sub- 
Eddington flows onto a stellar remnant black hole, but are 
not associated with galactic nuclei so are not powered by 
accretion onto a central supermassive black hole. If the Ed- 
dington limit holds, then the central object must be inter- 
mediate in mass between the population of < 20Mq stel- 
lar remnant black holes known in our Galaxy and the su- 
permassive black holes which form in the centres of galaxy 
bulges, and are therefore termed intermediate-mass black 
holes (IMBHs). However, there are many serious problems 
with such an interpretation, not least the in situ formation 
of the large numbers of such objects required to explain the 
populations of ULXs in star forming regions, that would re- 
sult in a n unfeasibly high fraction of mass being held in the 
IMBHs iKind 120041 ). Thus the majority of ULXs are most 
probably stellar remnant black holefl, accreting at super- 
Eddington rates (e.g. King et al. 2001; Watarai, Mizuno & 
Mineshige 2001; Begelman 2002; see also Roberts 2007 and 
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^ It is possible that stellar remnant black holes may range up to 
~ 90Mp, in mass if formed in low-metallicity environments, see 
IZampieri fc RobertsI bOOgll . 



references therein), though in the absence of a direct dy- 
namical mass measurement we cannot yet confirm this for 
any individual ULX. Possibly the best current IMBH can- 
didates amongst ULXs are those with the highest observed 
luminosities, and those with mHz QPO detections (although 
see below on the latter point). However, even amongst the 
brightest ULXs (the hyperluminous X-ray sources, or HLXs, 
e.g. Gao et al. 2003) there may be explanations other than 
an IMBH - for example, the HLX reported by Miniutti et al. 
(2006) has subsequently been revealed to be a background 
QSO (G. Miniutti, private communication), and the record- 
breaking HLX reported by Farrell et al. (2009) may turn 
out to be the stripped nucleus of a dwarf galaxy (Soria et 
al. 2009). 

The bulk of the ULX population therefore probably 
probes a new mode of accretion onto stellar-remnant black 
holes, different to the sub-Eddington states that are well 
studied from observations of black hole binaries (BHBs) in 
our own Galaxy (see e.g. the reviews by Remillard & Mc- 
Clintock 2006 hereafter RM06; Done, Gierlinski & Kubota 
2007 hereafter DGK07). Thus they can constrain models of 
super-Eddington accretion flows, which have wider impact 
than just ULXs. In the local Universe, these flows power 
objects such as the unique BHB GRS 1915-1-105, and the 
puzzling narrow line Seyfert Is. They also have cosmologi- 
cal signiflcance as super-Eddington flows are required in the 
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early Universe in order to quickly build up black hole mass 
to form the first quasars. ULXs can then give a template for 
understanding super-Eddington accretion flows, in a similar 
way to that of the BHBs for sub-Eddington flows. 

For BHBs the extensive, high quality monitoring data 
available means that their spectral variability as a function 
of luminosity can be studied for individual objects as well 
as for the class as a whole. This has given rise to a well 
defined picture, where at low luminosities the spectra are 
dominated by a hard F < 1.8 power-law tail extending out 
to ~ 100 keV. The thermal emission from an optically thick 
disc starts to become more important, increasing in both 
temperature and luminosity as the source brightens. This 
is correlated with a softening of the power-law tail, which 
pivots towards F ~ 2.1 (low/hard state: LHS). The spec- 
trum then makes an abrupt transition to the high/soft state 
(IfSS) where the disc dominates the spectrum, though this is 
generally accompanied by a weak non-thermal tail to high 
energies with F > 2.2. In this state the disc varies with 
L oc r** as expected for a constant inner disc radius. How- 
ever, at these high luminosities the non-thermal tail can also 
become strong, carrying a large fraction of the power, at 
which point the disc spectrum is strongly distorted by this 
Comptonisation (very high or steep power-law state: VHS) 
see e.g. the reviews by RM06; DGK07. 

By contrast, the much lower flux from ULXs means in- 
stead that there are typically only snapshot observations 
of individual objects, so putting together a picture of the 
spectral variability is much more difficult. Nonetheless, this 
has been attempted, though the limited quality of the data 
mean that often the spectra can be equally well described by 
power-law or disc models, which obviously confuses the spec- 
tral state identification (e.g . Glad stone & Roberts 2009). 
iKubota. Done fc Makishimal (|2002l ') used two ASCA obser- 
vations of IC 342 to track the spectral change of the two 
ULXs in this field. Both ULXs show a transition between a 
clearly curved spectrum which can be fit with a disc model, 
to one which can be fit by a power-law. Previous work has 
interpreted such changes as being from the HSS to LHS, 
and as this transition generally occurs at sub-Eddington lu- 
minosities then this would favour an IMBH. Kubota et al. 
(2002) instead showed that this could be interpreted as the 
much higher luminosity HSS to VHS transition, consistent 
with a stellar remnant black hole. 'Spectral transitions' from 
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ULXs have spectra which become harder as they become 
bri ghter, while others sho w the opposite trend. Recent work 
bv lFeng fc Kaared (|2009l ') and iKaiava fc PoutanenI (|2009l ') 
has begun to systematically analyse XMM-Newton data for 
individual ULX variability patterns, and interpret it in a 
physical framework. 

Another approach is to use higher signal-to-noise data 
from many different ULX to try to delineate the spectral 
changes as a function of luminosity, though this is compli- 
cated by the fact that ULXs are likely to be a heteroge- 
neous population of X-ray sources, in terms of black hole 
mass and viewing angle. Generally these spectra are fit with 
disc plus power-law models, and the derived disc temper- 
atures are sometimes low, ~ 0.2 keV, which (if taken at 



face value) implies a high black hole mass of ~ 10"^ Mq 
(e.g. Miller et al. 2003; Kaaret et al. 2003; Miller, Fabian & 
Miller 2004). While this initially supports the IMBH LHS 
interpretation, these higher quality data also clearly show 
that the accompanying power-law tail curves at the high- 
est energies, with a deficit of photons above 5 keV (Stob- 
bart, Roberts & Wilms 2006; Miyawaki et al. 2009). Such 
curvature is never seen in the BHB LHS at these low ener- 
gies, arguing against a simple IMBH interpretation for these 
objects on the basis of their X-ray spectra. On a similar 
theme, the few tens of mHz QPOs seen in some ULXs imply 
IMBHs if the observed QPOs are the direct analogues of the 
strong low frequency type-C QPOs seen in BHBs at 1-10 Hz 
(e.g. Strohmayer & Mushotzky 2009). However, we note that 
the likely super-Eddington, ~ 14M0 BHB GRS 1915+105 
can show QPOs at similarly very low frequencies (Morgan, 
Remillard & Greiner 1997, Hoel, Vaughan & Roberts 2009). 
Such QPOs are then clearly not an unambiguous signature 
of IMBHs. 

The highest quality X-ray spectra are actually much 
better fitted by disc plus low temperature, optically thick 
Comptonisation models (Stobbart et al. 2006; Roberts 2007; 
Gladstone, Roberts fc Done 2009). Using this description, 
[Gladstone et all (|2009l ') sorted the 12 best available XMM- 
Newton ULX spectra into a possible luminosity sequence, 
starting from a disc dominated, bright high/soft state at 
around the Eddington limit. For higher mass accretion rates 
they suggested an increasing fraction of the power is dissi- 
pated in a corona which covers the inner disc. This corona 
increases in optical depth and decreases in temperature as 
the source brightens, perhaps due to the increase in mass 
loss from the inner disc for these super-Eddington flows. At 
the highest luminosities the corona is strongly outflowing, 
giving rise to an increasingly large photosphere around the 
source which thermalises the energy to progressively lower 
temperatures (e.g. Kawashima et al. 2009; Begelman, King 
fc Pringle 2006; Poutanen et al. 2007). 

While plausible, such models must be tested on in- 
dividual objects where the mass and inclination of the 
black hole are fixed since ULXs seem likely to be a het- 
erogeneous population, as stressed above. SWIFT has been 
conducting a monitoring program for some ULXs in re- 
cent years (|Moore fc Milleij|2008l : iKaaret fc Fenell2009l '). In 
our present study, we use the data obtained from the 
most intensive of these, that of Holmberg IX X-1 (also 
historically known as M81 X-9). This source is located 
close to the dwarf galaxy Holmberg IX {d — 3.4 Mpc) 
and is the brightest X-r ay source within 20 a rcminutes 
of the nucleus of M81 (|La Parola et al.l [2toJ). This is 
also one of the ULXs that has been most widel y stud- 
ied in the past ten years (|La Parola et al.l 200 ll: IWand 
2005; 'Mil ler et al.1 12 004': 'Dewangan. Griffi ths fc Racl l2006l : 



Tsun oda et al.1 12008 : Kaaret fc Feng 2009|). However, all 
these either used hardness ratios to follow the spectral vari- 
ability, or at best fit disc plus power-law models. Here we 
choose to follow the spectral variability of Holmberg IX X- 
1 from the Swift monitoring data, but model this in de- 
tail using the disc and low temperature, optically thick 
Comptonisation which is required to f it the highest quality 
XMM-Newton data fr om this source (jStobbart et al.ll200a : 
[Gladstone et al.ll2009D . Section 2 describes how we extract 
and bin the intensity sorted spectra, while Section 3 gives 
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Figure 1. The Swift light curve and hardness ratio for Holmbcrg IX X-1. We show the 0.3 - 10 keV count rate as a function of observation 
time (top panel), and the corresponding hardness ratio (HR, as defined in the text; bottom panel). 



details of the spectral fitti ng results. Section 4 de scribes how 
these broadly fit into the [Gladstone et al.l (|2009l ) picture of 
spectral variability described above, but also shows the need 
for better data to understand these trends. 



2 DATA 

We use only those Swift XRT observations that are pointed 
at Holmberg IX X-1, so we can directly compare instrument 
count rates rather than having to correct for the different 
instrument response of off-axis pointings. This gives a total 
of 75 pointings, each typically with ~ 1700 s exposure, from 
Swift pr ogram 00090008 , cover ing the whole monitoring pro- 
gram o f lKaaret fc Fend (|2009l ) up to July 2009. In addition, 
we also use 7 observations from Swift program 00035335 
during 2006 - 2007. These are longer observations (up to 
20ks) but have a much sparser coverage of the light curve. 

We extract the photon counting mode of the Swift 
pointed observation data using the Suii/f-specific ftool 
xrtpipeline to obtain clean event files. We use standard 
data selection criteria, with minimum elevation angle (ELV) 
set to be larger than 45°, the bright Earth (BR_EARTH) 
angle set to be above 120°, minimum Sun and Moon angle 
(SUN_ANGLE and MOON_ANGLE) above 45° and 30°, re- 
spectively, and pointing direction within 0.08° of the source. 
We then extract spectra for both source and background 



from grades 0-12 within a 47 arcsec (20 pixels) radius 
region centred on the ULX. A background spectrum was 
extracted from a same-sized aperture on the edge of the 
field-of-view. 

Fig. 1 (upper panel) shows the resulting light curve for 
these data over the 0.3 - 10 keV band. Plainly the source is 
persistent, but variable by a factor of 3 - 4 in intensity. We 
also show the hardness ratio (HR) for these data, defined 
using the standard Swift XRT hard and soft energy bands 
as the rate from 1.5 - 10 keV divided by the rate from 0.3 
- 1.5 keV. The lower panel on Fig. 1 shows that this HR 
does vary, but with large uncertainties associated with the 
short observations. We show the variation of HR with inten- 
sity in Fig. 2, where the 7 data points from the 2006 - 2007 
observations are shown in (magenta) filled stars while those 
from the monitoring program are shown as (cyan) crosses. 
This shows that the spectrum hardens slightly as the lu- 
minosity increa s es, as seen in the analysis of these data by 
iKaaret fc Fend (|2009l ) (although their hardness ratio is de- 
fined slightly differently, as a ratio of 1 - 8 keV/0.3 - 8 keV). 

We co-add observation data in order to increase signal- 
to-noise. We define 5 count rate ranges as detailed in Table 
1 to have roughly (within a factor 3) equal numbers of total 
counts as shown by the (red) dotted lines on Fig. 2. Data 
within these ranges are consistent with the same HR except 
for SWIFT 5, the highest intensity bin, where the long ob- 
servation in 2006 has a slightly but significantly softer spec- 
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Figure 2. The hardness ratio-intensity diagram for Holmborg 
IX X-1. We show each individual observation data point from 
the monitoring programme as a cyan cross, while the data points 
from the earlier, longer observations are indicated by magenta 
filled stars. The count rate ranges used for co-adding the data are 
delineated by the red dashed vertical lines, and the values from 
these co-added data are shown by black circles. 



trum than the ones from the monitoring campaign. These 
data cannot be co-added as their spectra are not consistent. 
Hence we include the 2006 data only for this intensity bin 
as this has higher statistics. 

We add the spectra of each count-rate group using the 
addspec FTOOL. As recommended, the response files are 
weighted with their counts. The differences between the aux- 
iliary response files (*.airf) over all data sets in each count- 
rate range is less than 10%. We include the hardness ratios 
from these co-added spectra on the hardness ratio-intensity 
plot (Fig. 2) as the large open circle (black) points. For the 
purposes of spectral fitting, the data were binned with at 
least 20 counts in each bin. 

We also include the long exposure (114ks) XMM- 
Newto n data of Holmberg IX X-1 from [Gladstone et al.l 
l|2009l ). hereafter termed XMMl (observation ID 
0200980101, taken on 2004-09-26). The spectrum of 
this XMM- Newton data is r educe d by using the same 
method as in [Gladstone et all (|2009l ). These data have very 
similar 0.3 - 10 keV flux to that of the lowest luminosity 
Swift spectrum. However, its hardness ratio (when calcu- 
lated independently from the instrument response using 
model fluxes in the bands) is signiflcantly softer, closest to 
that from SWIFT 4. This clearly shows that the variability 
is not simply a function of luminosity, and that there is 
extra complexity present. 

There are two other XMM-Newton data sets available 
in the archive (hereafter XMM2 and XMM3). These are all 
between XMMl and SWIFT 2 in intensity, so do not span 
much of a range in flux. They also span the range in shape 
between XMMl to SWIFT 2, so do not add to the spectral 
variability but we include them for completeness. 



3 MODEL FITTING AND RESULTS 

In some earlier studies, XMMl data was conventionally anal- 
ysed using the 'standard' disc plus power-law model (e.g. De- 
wangan et al. 2006; Winter, Mushotzky & Reynolds 2007). 
While this conventional model provides a reasonably accept- 
able fit, the spectrum curves at high energy such that a 
power-law tail overpre dicts the observations above 5 keV 
({Gladstone et al.1l2009l l. Instead, a disc plus low tempera- 
ture (and hence optically thick) Comptonisation spectrum 
fits this spectral shape, as the low electron temperature gives 
the low energy rollover in the Com ptonised flux. Hence we 
use this model to flt all our spectra. [Gladstone et al.l (|2009l ) 
used some very sophisticated disc and Comptonisation mod- 
elling to fit their sample containing the longest exposure 
XMM-Newton ULX data available in the archive. Here we 
have more limited signal to noise in the Swift spectra so we 
choose the simplest possible model for disc plus low temper- 
ature Comptonisation, namely a diskbb plus comptt model, 
and set the seed photon temperature equal to the disc tem- 
perature. 

We assume a HI column density of 4.06x10^° cm~^ 
in th e direction of Holmberg IX X-1 (|Dickev fc LockmanI 
and we allow intrinsic absorption in the host galaxy 
to be a free parameter. We fit our 8 spectra (3 XMM-Newton 
plus 5 Swift spectra) simultaneously in xspec, allowing all 
parameters to be free except for the extra-galactic absorp- 
tion column which is tied across all 8 datasets. This gives 
a reasonable fit, with xi ~ 3512/3461 for a best fit value 
of Nh = (1.04 ± 0.02) X 10^^ cm~^ The derived disc and 
Comptonisation parameters are given in Table 2, where the 
fluxes are the unabsorbed (intrinsic) flux from the diskbb 
plus comptt model from 0.3 - 10 keV. The error of the flux is 
calculated by using flux err command in xspec. The bolo- 
metric luminosity is calculated from the unabsorbed flux in 
0.01 - 100 keV range, assuming disc geometry. 

All the spectra are well fitted by a low temperature 
Comptonising corona, that is optically thick. However, the 
electron temperature only has a firm upper limit in XMMl, 
XMM2, SWIFT 3 and SWIFT 4, so in the other datasets 
the rollover is not significantly detected. Nonetheless, the 
fact that some of the data (especially XMMl, see Fig. 5 of 
Gladstone et al. 2009) do require this rollover shows that 
we should use this model for all the data for consistency. 
Similarly, the moderate signal-to-noise of the SWIFT data 
means that the disc component is not strongly required in 
any of them (Ax'^ ^ 7 for one additional degree of freedom) 
but we again include it for consistency. 

Fig. 3, left panel, shows the derived spectra for XMMl 
and the SWIFT spectra, after correcting for absorption in 
both our Galaxy and the host galaxy, and deconvolving from 
the instrument response. XMM2 and XMM3 are not in- 
cluded on this plot as they overlap with XMMl, SWIFT 1 
and SWIFT2. Instead, the three XMM spectra are shown 
separately in the right panel. 

These plots clearly show the complex nature of the spec- 
tral variability. XMMl has a marked inflection at 1 - 2 keV, 
separating the emission into two peaks which are flt by the 
cool disc and low temperature Comptonisation at low and 
high energies, respectively. XMM2 and SWIFT 2 have a sim- 
ilar but less pronounced feature. None of the other spectra 
appear to have two spectral peaks. This is especially in- 
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Table 1. The table details the datasets we use from Swift to make the co-added spectra with higher signal- to- noise and 
XMM-Newton. The prefix M and L on Swift denotes 000900080 and 000353350, respectively. 



Name count rate total counts Swift and XMM-Newton Obsid 



SWIFT 1 


0.00-0.15 


4882 


M04, M12, M13, M14, M17, M18, M26, M28, M34, M45, M46, M47, M65, M66, 
M67, M68, M69, M70, M71, M72, M73, M75, M76, M77, M78, M79, M81, M82. 


SWIFT 2 


0.15-0.20 


11733 


Mil, M15, M19, M22, M25, M29, M30, M31, M32, M33, M36, M38, M41, M44, 
M48, M49, M50, M52, M54, M55, M56, M57, M58, M64, LOl, L02, L07, L09. 


SWIFT 3 


0.20-0.25 


7143 


M02, M03, M05, M07, M20, M21, M23, M24, M39, M43, M59, M60, M61, M62, 
M63. 


SWIFT 4 


0.25-0.30 


6551 


M08, M09, MIO, M40, M42, L05, L06. 


SWIFT 5 


0.30-0.35 


6773 


L04 


XMMl 


1.53 ±0.005 


114000 


0200980101 


XMM2 


1.80 ±0.02 


12800 


0112521001 


XMM3 


2.07 ±0.02 


16000 


0112521101 




Figure 3. Spectral data and corresponding best fitting model. The data arc de-absorbed and unfolded from the detector response. Left 
panel from bottom to top: XMMl (light blue), SWIFT 1 (magenta), SWIFT 2 (red), SWIFT 3 (green), SWIFT 4 (blue), and SWIFT 
5 (cyan), respectively. Right panel from bottom to top: XMMl (light blue), XMM2 (orange), XMM3 (light green). The SWIFT data 
groups are plotted to have at least 8cr significance or arc grouped in sets of 10 bins, while the XMM-Newton data is plotted to a minimum 
20(T significance or 20 bins for clarity. 



triguing for SWIFT 1, as this has very similar overall flux 
to XMMl yet is plainly different in shape, with more flux 
between 1-2 keV filling in the infiection, and less flux below 
0.7 keV. This difference could instead be due to increased 
absorption in SWIFT 1 rather than a real spectral change for 
the constant absorption column assumed here. However, the 
low energy spectral shape of SWIFT 1 is similar to that of 
the other SWIFT spectra, so it seems more likely that there 
is a real spectral change between the SWIFT spectra and 
XMMl. This is not from the different instrument responses 
as XMM3 has a similar low energy slope to SWIFT 2. Thus 
it seems most likely that there are real spectral changes for 
little or no change in luminosity. 

At higher energies there is more variability. The high en- 
ergy rollover clearly indicates that the Comptonising corona 
has a lower electron temperature in SWIFT 3 and 4 than in 
the lower luminosity spectra (XMMl - 3 and SWIFT 1 and 
2), but this trend is not continued to the highest luminos- 



ity data (SWIFT 5), where the rollover again moves up in 
energy to close to the edge of the bandpass. 

We plot the model parameters in Fig. 4 as a func- 
tion of the unabsorbed, bolometric luminosity. The XMM- 
Newton data are plotted as filled stars: XMMl, XMM2, 
and XMM3, from left to right, respectively. The disc is 
consistent with having the same temperature in all spectra 
ixi ~ 3512/3461). Alternatively, it is consistent with having 
the same normalisation in all spectra, but tying both tem- 
perature and normalisation give a significantly worse fit at 
X? ~ 3562/3475. This gives the almost constant low energy 
flux described above. 

The electron temperature is clearly changing, decreas- 
ing with increasing luminosity through XMMl - 2 and 
SWIFT 1-4, except for XMM3 where we cannot have a 
good constraint on the upper limit of the electron temper- 
ature, even worse than SWIFT 5 (and thus its bolomet- 
ric flux has likely been overestimated in Fig. 4 and Table 
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Table 2. Spectral fitting results. The flux is in units of 10~^^ erg cra~^ while the bolometric 
luminosity, Lboh is in 10^^ erg s~^. See text for the range of count rate for each SWIFT data group. 
The errors quoted here and elsewhere in the paper arc the 90% confidence interval for one interesting 
parameter. 



Group 




Norm.^/^ 




T 


Flux_^ 


b 




SWIFT 1 


0.203 M 

— u.uoy 


5.04+1-^^ 

— 2.0Z 


3.01+l\l 

— i. ii 


6.86+1-11 


6.22+0-^!» 

— 4.4 


6.63 


189.6/171 


SWIFT 2 






q 1 1 +58.7 
•-'•-^-^-0.828 


6-6015.04 


8.70t°:«f 


9.27 


337.8/327 


SWIFT 3 


r, r,ic:+0.105 
'^■^^''-0.059 




±.uo_o 223 


^■^"-1.17 


9.36+°:?f 


9.41 


226.1/245 


SWIFT 4 


„ r,c:o+0.359 

U.ZOZ_(, ggg 




-, .7+0.221 
-0.161 




19 0+0. 600 
i^-»-1.80 


11.1 


232.1/236 


SWIFT 5 


0.488^0;^? 


9 Qcr+3.85 
^•'-'•^-0.964 


9 47+38.4 
-0.813 


8 1 5+42-2 
''■^^-6.32 


18.0lif, 


17.3 


227.1/244 


XMMl 


n 907+0.011 

'^■^'^'-0.05 


5 (57+0.434 
^•°'-0.395 


9 44+0.158 
^■**-0.140 


8-76tg;g^ 


rr 05+0.270 
°-^>'-0.570 


6.45 


1349.3/1286 


XMM2 




4.98tllt 


9 01 +0.901 
^"^^-0.309 


8-0911:82 


fj fi5+0.460 


6.92 


435.7/447 


XMM3 


159+0.056 




^■°'-6.61 




7 Qi +0.890 
' •''^-6.23 


11.1 


514.2/506 



Unabsorbed flux in the fitting energy range (0.3 - 10 keV). 
^ The bolometric luminosity is calculated from unabsorbed flux in 0.01 - 100 keV range. 



2). The similar spectral slope of all the data in the 2 - 
4 keV range dominated by Comptonisation requires a simi- 
lar Compton y parameter, defined as y = At^O — constant, 
where 6 = kTcjra^c? ^ so the lower temperature implies an 
increased optical depth. A lower electron temperature also 
brings the rollover further into the observed bandpass, so 
it is better constrained. Hence the optical depth (derived 
from the observable parameters of spectral slope and elec- 
tron temperature) has correspondingly smaller errors also. 
However, the highest luminosity spectrum does not follow 
this trend. The electron temperature increases again, so the 
optical depth drops. 



4 DISCUSSION 

[Gladstone et all (|2009l ) suggested that the different spec- 
tral shapes seen in their sample of ULXs could be set into 
a sequence such that the corona became increasingly op- 
tically thick and lower temperature, extending over more 
of the inner disc, as the mass accretion rate M increased 
beyond Eddington (MEdd)- Plausibly these changes could 
be caused by the increasingly strong winds which are ex- 
pected to be driven by such super-Eddington flows. At even 
higher M /MEdd the source becomes embedded in an ex- 
panding photosphere, l owering the temperat ure still further 
iBegelnian et al.ll2006l : iPoutanen et aLlboO?! ). 

Our results show some support for this picture. Figs. 3 
and 4 show that the electron temperature and optical depth 
in the corona generally decrease with increasing luminosity 
except for the brightest spectrum. However, they also high- 
light the fact that the spectral variability does not uniquely 
track the observed luminosity. SWIFT 1 and XMMl have 
very similar fluxes, but their spectra are subtly different. 
Such lack of one-to-one correspondence between spectral 



S o 



Z 





1 ' 


1 1 1 1 1 1 


1 1 E 






H ^ 

H — H 


— 1 1 1 1 1 1- 




1 r 

* 

'- 1 


-1 1 hi 


+ -1 


1 1 1 1 1 1 H 


1 ' ' 

1 1 '- 


'_ i 


-1 1 1- 

t 


H ^ 

1 \ 

1 


1 1 1 1 1 1- 







10 15 

L.^.ClO-^ergs s-') 



20 



Figure 4. Dependence of the best fitting parameters on the bolo- 
metric luminosity. From top to bottom: disc temperature kT\^ 
(keV); square root of diskbb normalization (which is oc disc inner 
radius); electron temperature in the corona kTc (keV); and the 
electron scattering optical depth r. Filled stars are XMM-Newton 
data: XMMl, XMM2, and XMM3 from left to right, respectively. 
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shape and luminosity is a well-known feature of BHBs for 
the transition between the low/hard and high/soft states. 
However, this hysteresis in BHBs is most probably con- 
nected to non-equilibrium effects caused by the rapid rise 
in mass accretion rate in the transient out bursts of these 
systems ([Gladstone. Done fc Gierlinskill2007l ) , whereas most 
ULX are persistent sources. BHBs can also however show an- 
other type of non-unique relation of spectral shape to lumi- 
nosity, where the strongly Comptonised very high state can 
be see n at the same luminosity as a standard disc dominated 
state jGierlinski fc Don3l2004l . and the references therein). 
It is not at all clear what determines which of these high 
mass accretion rate states is shown at any time. The dif- 
ference between SWIFT 1 and XMMl may then be con- 
nected to the same (unknown) mechanism. Alternatively it 
may be connected to the additional complexity expected 
from super-Eddington flows. Numerical simulations show 
that such flows are highly dynamical (Ohsuga 2006; 2007; 
2009; Kawashima et al. 2009; Takeuchi, Mineshige & Ohsuga 
2009). Hence changes in the disc/wind geometry along the 
line of sight could result in complex changes in the ob- 
served flux which do not reflect the intrinsic luminosity of 
the source. However, the complex structure in accretion flow 
shown by the numerical simulations mentioned above may 
have been smoothed out within the observation timescale of 
our present study (of the order of days). Alternatively, we 
may consider an intermittent outflow (e.g. Kato, Mineshige, 
& Shibata 2004). 

It is interesting to compare the variability behaviour of 
Holmberg IX X-1 to other ULXs modelled wi th similar cool 
disc p lus Comptonisation models. Specificallv. lRoberts et al] 
(|2006l ) explored this for the ULX NGC 5204 X-1 based on 
a C handra monitoring pro gramme and XMM-Newton data, 
an d iFeng fc Kaaretl (|2009l ') use this model to compare XMM- 
Newton data for several ULXs over all available epochs with 
sufficient data quality. Interestingly, these works show that 
IC 342 X-1 behaves similarly to Holmberg IX X-1; its coronal 
electron temperature appears to drop as it gains in luminos- 
ity. In contrast, both NGC 5204 X-1 and Holmberg II X-1 
behave differently. In these ULXs the electron temperature 
increases with increasing luminosity, with a corresponding 
drop in optical depth. Perhaps most interestingly, there ap- 
pears to be a correlation between these spectral changes and 
the position of the sources in the spectral sequence for ULXs 
shown as Fig. 8 of Gladstone et al. (2009). Both Holmberg 
IX X-1 and IC 342 X-1 appear in the central portion of the 
sequence, where the spectra appear with a strong inflection 
and the energy density peaks in the higher energy, optically 
thick, Comptonised spectral component. Holmberg II X-1 
and NGC 5204 X-1, on the other hand, appear further along 
the sequence, where the energy density peaks in the softer 
part of the spectrum and Gladstone et al. (2009) hypothesize 
that the spectra are becoming photosphere-dominated due 
to the increasingly important effects of a strong wind. This 
apparent correlation between spectral shape and variabil- 
ity data supports the conclusion of Gladstone et al. (2009) 
that different physical processes underly the two distinct 
spectral shapes. However, we emphasize that this tentative 
result is concluded from the current, sparse data available 
for ULX spectral variability studies; clearly new dedicated 
observational programmes are required to investigate these 
phenomena in more detail. 



5 CONCLUSIONS 

We use the multiple Swift observations of Holmberg IX X-1, 
co-added in different intensity ranges, together with a long 
XMM-Newton observation of this source, to make the most 
sensitive study to date of spectral variability in any individ- 
ual ULX. We show that the variability is subtle and rather 
complex, such that it is not well characterised by a single 
hardness ratio. Instead, we do full spectral fitting with a 
disc plus thermal Comptonisation model. There is a general 
trend in Holmberg IX X-1 for the Comptonising corona to 
decrease in temperature and increase in optical depth as the 
so urce luminosity i ncreas es. This supports the suggestion 
of [Gladstone et al.l (|2009l ') that an increase in mass loss in 
winds from the inner disc as the source becomes more super- 
Eddington leads to increased material in the corona, thereby 
increasing the optical depth. The acceleration mechanism is 
then required to share the coronal power across more par- 
ticles, so the temperature drops. We also note that there 
appears to be a tentative correlation between the observed 
changes in the coronal parameters with luminosity, and the 
type of spectrum displayed, over the small number of ULXs 
with sufficient quality data to begin to explore spectral vari- 
ability. 

However, the data for Holmberg IX X-1 also show more 
complex behaviour. It is clear that the spectral shape is 
not uniquely determined by the observed luminosity. At the 
same luminosity the spectra can look subtly different, and at 
different luminosities the spectral shape can be the same. It 
may indicate that some fraction of the observed variability is 
due to time dependent line of sight column changes rather 
than intrinsic changes. Alternatively, this may be related 
to the equally puzzling lack of one-to-one correspondence 
between spectral shape and luminosity seen in the BHBs at 
high (but apparently sub-Eddington) mass accretion rates 
where the source can show either high/soft or very high state 
at the same luminosity. Whatever the reason, this severely 
complicates any attempt to understand the origin of the 
spectral variability. It is therefore very evident that if we are 
to make progress in understanding ULX spectral variability, 
and the physical processes that underly these behaviours, 
then the acquisition of new observational data composed of 
good quality ULX spectra over a range of timescales and 
source fluxes is absolutely imperative. 
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